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Abstract 

It is known that similar physical systems can reveal two quite different ways of behavior, either coarsening, which creates a 
uniform state or a large-scale structure, or formation of ordered or disordered patterns, which are never homogenized. We present 
a description of coarsening using simple basic models, the Allen-Cahn equation and the Cahn-Hilliard equation, and discuss the 
factors that may slow down and arrest the process of coarsening. Among them are pinning of domain walls on inhomogeneities, 
oscillatory tails of domain walls, nonlocal interactions, and others. Coarsening of pattern domains is also discussed. 
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1. Introduction 

For many decades, the nonlinear development of instabilities in physical systems was an object of extensive in¬ 
vestigations. The most spectacular consMuences of instabilities are the appearance of ordered spatially non-uniform 
structures (pattern formation, see 111] - lU) or irregular motions (spatio-temporal chaos flj) under uniform external 
conditions. However, there is one more scenario of an instability development; that is a gradual growth of the char¬ 
acteristic scale with time (coarsening) II, 0]. Different evolution scenarios can take place in rather similar physical 
systems. 

As an example, let us consider the phase separation in binary alloys that consist of two kinds of atoms, A and 
B, with volume fractions <^a(x, f) and 0 b(x, t), respectively. There exists a temperature Tc such that for T > Tc the 
components are mixed, i.e., the order parameter (p(x, t) - (pA('t^, 0 - t) vanishes everywhere, while for T < Tc 
they are separated, i.e., there exist two thermodynamically stable phases, one with <p > 0 (“A-rich phase”) and the 
other with f < 0 (“B-rich phase”). A mathematical model of that phenomenon has been suggested by Cahn and 
Hilliard I?!]. Under the simplest assumption of a constant mobility, the kinetics of the phase separation is described 
by the following non-dimensional equation (the Cahn-Hilliard equation), 

+ ( 1 ) 

The uniform phase <^ = 0 is unstable while the uniform phases 0 = +1 are stable. The instability of the phase f - 0 
creates a mosaic of islands of both stable phases. The size of these islands (domains) grows due to coarsening, which 
eventually leads to a complete separation of stable phases S, 0. 

A diblock copolymer, which consists of monomers A and B with reduced equal local densities (pA and fs, is quite 
similar to a binary alloy. The basic difference is the existence of a long-range interaction of monomers 0 - in 
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which provides an additional term in the evolution equation for the order parameter: 

(pt - v^(-0 + - v^<^) - r^, r > 0, {(p) - 0. 


( 2 ) 


There are no other spatially uniform stationary solutions except (p - 0. Therefore, when the latter solution is unstable 
(at T < 1/4), a transition to a non-uniform state is unavoidable 111 ill . At small T, the initial evolution of disturbances 
is similar to that in the Cahn-Hilliard equation, but it is stopped when stripes with a definite pattern wavelength are 
created 1121. 

The existence of long-wave linear instability and multiple homogeneous states does not guarantee the creation of 
spatially uniform domains through coarsening. As an example, let us discuss the nonlinear dynamics governed by the 
one-dimensional Kuramoto-Sivashinsky equation. 


(pt - -<pxx - (pxxxx - (<P^)x, {(p) - 0, 


(3) 


which is used for the description of instabilities in reaction-diffusion systems 111311 . 111411 . instabilities of flame fronts 
M, and film flow instabilities ifT^ . In that case, the linear part of the equation is identical to that of the Cahn-Hilliard 
equation, and any constant solution, <p - tpQ, satisfies equation (l3]l. However, all solutions corresponding to uniform 
states are unstable. The Kuramoto-Sivashinsky equation is a paradigmatic model of the spatio-temporal chaos iQ; 
stable periodic patterns are also possible ES, 10, but the attraction domain of that regime is small. Generally, the 
way of the instability development depends significantly on the details of the system nonlinearity and symmetry lIlSll . 

ill. 


Note that pattern formation and coarsening are not incompatible phenomena. Let us return to model (|2]i that de¬ 
scribes formation of stripes. Because of the rotational invariance of the problem, the orientation of stripes is arbitrary. 
Initially, a disordered system of stripes is developed from random initial conditions, and then the mean size of ordered 
domains grows with time, i.e., domain coarsening takes place for differently oriented stripe patterns rather than for 
different uniform phases 0, il. 

One can see that the interplay between coarsening and pattern formation is nontrivial, and it is the subject of the 
present chapter. Let us emphasize that here we discuss only dynamic models of coarsening and pattern formation, 
which do not include any kind of noise. The phenomena caused by thermal fluctuations are considered in other papers 
of the present issue and in the comprehensive book 0, where the reader can find many additional references on that 
subject. 


2. Coarsening in one dimension: dynamics of domain walls 


When considering the kinetics of coarsening, one has to take into account the following basic factors. 

1. The existence of a Lyapunov functional. If the temporal evolution of an n-component order parameter </>,, 
i - 1..... n is described by a gradient evolution equation 


k=l 




(4) 


where F - J L(tp, tp^, ■ ■ )dx is the Lyapunov functional of the system, and Dij is a positive definite matrix, then 


F,^ 




6F 6F , 

Dill - dx < 0. 

6(pi 5(pk 


(5) 


In the case of an equilibrium phase transition, the existence of the Lyapunov functional (free energy) is the conse¬ 
quence of the thermodynamics. The nonlinear development of instabilities in systems far from equilibrium in some 
cases is also governed, at least approximately, by potential systems of equations that possess Lyapunov functionals. 

2. The existence of a conservation law. In the case where the order parameter is a density of a conserved quantity 
(e.g., the number of molecules), the evolution equations have a divergence form. 


(pi,t - -V ■ J/, i- 1,..., A, (6) 

where the flux J, is a function of the order parameter and its derivatives. 

In order to understand how both factors influence the coarsening kinetics, let us consider a number of examples. 
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2.1. Non-conserved order parameter 


2.1.1. Allen-Cahn equation 


We start with the Allen-Cahn equation 1122 1 


(jit - 4>xx + (j) - <f)^. 


(7) 


which describes a phase transition in the absence of a conservation law for the order parameter (p. The physical 
interpretation of that model can be as follows [@]: (^ is the spontaneous magnetization directed along the definite axis 
(due to the crystal anisotropy). Also, the Allen-Cahn equation is the simplest example of an order parameter equation 
for a pattern forming system far from the thermodynamic equilibrium. Let us consider the onset of convection in 
a horizontal cylinder heated from below. Above the instability threshold, the rotation of the liquid in the transverse 
section of the cylinder can be either counterclockwise (0 > 0) or clockwise {(p < 0). The temporal evolution of the 
order parameter (p{x, t) {x is the coordinate along the axis of the cylinder) is governed by the Allen-Cahn equation 

Q. 

Equation O has a Lyapunov functional, 

/ \ I 

L{x, t)dx, L(x, t) - -<pI -I- -(0^ - 1)^ > 0. (8) 

Its derivative 


j 

rldL dL \ i 

^ y r 7 




{-(p + (p - <pxx)ftdx - - 1 dx(p, < 0 

(9) 


is non-positive, hence the Lyapunov functional decreases monotonically with time until a stationary state is reached. 
It is obvious that equation (|7]) has three fixed points corresponding to uniform phases, (p - 0 (unstable paramag¬ 
netic phase) and (p - ±\ (ferromagnetic phases with opposite orientations of the magnetization); the latter solutions 
correspond to the absolute minima of the Lyapunov functional, F -0. The separatrices 


4 >± 


- ± tanh 


V2 ’ 


^ = const 


( 10 ) 


{kink and antikink) describe domain walls separating semi-infinite domains with different signs of the order parameter. 
Note that 

(p± ~ +{1 - 2 exp[+(v - f) V2]) as V ^ +oo. (11) 

The contribution of a domain wall to the Lyapunov functional is Eq = 2 V2/3 > 0. 

Other stationary solutions of (|7]i describe spatially periodic patterns which can be considered as periodic arrays 
of domains with alternating signs of (p. They are expressed through the elliptic Jacobi function. 

Let us impose an infinitesimal disturbance 0(v, f) on the stationary solution (p{x) and consider its temporal evolution 
governed by the linearized problem, 

— ^xx + (1 ~ 30^(x))^, 1^1 < oo as X ^ + 00 . (12) 


For normal disturbances p(x) exp(crf) we obtain an eigenvalue problem, 

crp — H- (1 — 3(p^)p, |^(+oo| < oo. (13) 

It is obvious that solution 0 = 0 is unstable: cr(k) = 1 - P for 0 = and the solutions 0 = +1 are stable: 
cr(k) = -2 - tr for the same 0. The domain wall solutions (fTOl i are neutrally stable. One can show that all the 
spatially-periodic solutions mentioned above are unstable. 

Let us discuss now the temporal evolution of the system, when the initial state is the unstable phase <p - 0 with 
a certain initial, spatial disordered, perturbation. According to the stability analysis presented above, the final state 
should be uniform or contain a single defect, a domain wall. However, it is clear that the decomposition of the phase 
0 = 0 can produce numerous domains with alternating signs of (p. Such a state can be characterized by a certain 
density of defects decreasing with time, or by a mean domain length which grows with time. At late stages of the 
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evolution, when the typical distance between domain walls is large the analysis can be done by means of asymptotic 
methods 12^ . One considers a set of domain walls (doll of alternating signs centered at i — 1 , 2 ,..., 

and slowly moving because of their interaction. By means of asymptotic expansions, the original nonlinear problem 
is transformed into an infinite system of inhomogeneous linear equations. Their solvability conditions determine the 
following equations of motion for the centers of domain walls 12611 : 


2 V 2 . _ dU 


(14) 


where 

~ " -8 V2exp [- V2(^,- - . (15) 

i 

Thus, the domain walls attract to each other according to the exponential interaction law (fT4l l which reflects the 
exponential domain wall asymptotics (fTTT i. 

System (fT4l i has a family of stationary solutions 


^j-a + jl, 7 - 0 ,+ 1 ,+ 2 ,... (16) 

corresponding to periodic patterns with the spatial period 21. The attractive interaction makes all these solutions 
unstable. 

Consider the interaction of a pair of domain walls. Setting /(f) = ^ 2(0 - we find that the distance between 
domain walls is governed by the equation 

^/'(f) = -32e-'^. (17) 

If 1(0) = /o » 1, the solution is 

/ = /o+-^ln(l-48e-^'“f). (18) 

The distance between two domain walls becomes of 6/(1) at 

f~fo = ^e^'“+0(l). (19) 

Finally, the domain walls reach the distance of order 6/(1) and annihilate. During the time interval f, only the domain 
walls with the original separation greater than ~ ln(48f)/ V2 can survive. Thus, a logarithmic coarsening takes place. 
In a large but finite system with the length L, the number of domain walls N(t) < Lll(t). 

As an example of the situation when two locally stable phases are energetically non-equivalent, let us consider a 
system with the Lyapunov functional 


F[<f>{x)] 


f 


2 ^^ 4(1 




h(f> 


dx 


which corresponds to the dynamic equation 


( 20 ) 


+ (p - (fr + h. 


( 21 ) 


In the case of a magnetic system, the last term in the expression (l20l i describes the influence of an external magnetic 
field, which makes the orientation of the magnetization in the direction of the field preferable. The uniform stationary 
states satisfy the equation 

+h^0. (22) 

In the interval -/z, < h < /z,, equation (l2^ has three solutions: stable solution 4> - ‘P+ > V3, another stable 

solution 0 = < -1/ V 3 , and an intermediate unstable solution -1/ V3 < 0o < 1/ V3- For h > 0, the phase with 

(p - <p+ > 0 is stable, and the phase with (p - (p^is metastable. 


4 





/ Comptes Rendus Physique 00 (2015) 1 414\ 


5 


If h is small, domain walls are described by formulas (fTOl) at the leading order. By means of an asymptotic analysis, 
one can find that the motion of each domain wall is governed by the equation | 


= , 2 /. 

3 dt 


(23) 


(the upper sign is for a kink and the lower sign is for an antikink). This motion creates a coarsening process, which 
is significantly faster than that in the case of energetically equivalent phases. That process leads to the elimination of 
the metastable phase. 

2.1.2. Fractional Allen-Cahn equation 

The logarithmic law of the one-dimensional coarsening is caused by the exponentially weak interaction, due to the 
exponential asymptotics (fTTT) characteristic for solutions of partial differential equations. However, the basic equations 
governing the natural phenomena are often integro-differential equations rather than partial differential equations. For 
instance, the free energy density of a fluid depends on the fluid density in a nonlocal wa v 12711 . The conventional van 
der Waals’ expression for the fluid free energy which contains a squared density gradient 1281] corresponds to a certain 
asymptotic limit of the basic nonlocal expression. The asymptotics of a domain wall (i.e., the gas-liquid boundary) 
in the framework of local and nonlocal models are strongly different; while a local model predicts an exponential 
decay, a nonlocal model suggests a power-law decay 12^ . Fronts that cannot be governed by local partial differential 
equations have been found also in studies of transitions with long range interactions 1291] . vacancy diffusion and 
domain growth in binary alloys S, and ordering kinetic on fractal structures j^l]] ( for the latter subject, see ll^ l. 
As an example let us discuss the front propagation in systems with superdiffusion 1331] . While the normal diffusion is 
associated with a Gaussian non-correlated random walk of particles, the superdiffusion is observed in non-equilibrium 
systems with an algebraically deca ying jump length distribution, where the central limit theorem is no valid. Among 
the examples are wave turbulence il34] . transport in porous media (3, and forage trajectories of animals 1^ . One 
can use a superdiffusive generalization of the Allen-Cahn equation. 




, 1 <r<2. 


Here Dj. denotes the fractional Riesz derivative, which can be defined by its action in the Fourier space: 




(24) 


(25) 


where F is the symbol of the Fourier transform. Equation (l24l) has a Lyapunov functional 1371] . In a contradistinction to 
the exponential asymptotics (fTTT) of the domain walls in the local PDE (|7]i, the domain walls in the integro-differential 
equation (l24li has an algebraic tail. 


2r(2 - y) 

That leads to a power law for the kink-antikink attraction, 

i'(t) = -cn 

(cf. (HU), and for the temporal decay of the number of domain walls on a finite spatial interval, N ~ t OTT . 


(26) 


(27) 


2.1.3. Non-potential systems 

Eor systems far from equilibrium, e.g., in the case of longwave instabilities of flows, the Lyapunov functional 
generally does not exist. Nevertheless, coarsening may take place if the domain walls are described by monotonic 
functions. As an example, let us mention the amplitude equation that governs the fixed-flux convection in a tilted slot 

Q: 

<Pt — <Pxx + ^ + 2a(p(px, O' > 0. (28) 

Because the symmetry x ^ -x is violated, the kink and antikink domain walls have different widths: 

(^±(x) = tanhyS±(x - f), p± - ^(a ± -H 2). 

5 


( 29 ) 
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The interaction of domain walls is attractive but asymmetric. For a pair which consists of a kink with the center in the 
point (f) and an antikink with the center in the point ^2(0^ the equations of motion look as 

exp(-2|/?_|(^2 - ^i)), ^2 = F exp(-2yS+(f2 - ^i)), (30) 

where F+ 4^ F-. All the periodic stationary solutions (patterns) are unstable, and a logarithmically slow coarsening 
takes place, similarly to the case of the Allen-Cahn dynamics. 


2.2. Conserved order parameter 
2.2.1. Potential systems 

Let us return to the Cahn-Hilliard equation ([TJ, which can be written as 

(jit + jx = 0, = ((t>xx + <f>- cf>^)x', -1<X<1. (31) 

Being a generic nonlinear equation governing longwave instabilities in the presence of the conservation law ifl^ . that 
equation has been revealed in numerous problems of different physical nature, including secondary flows produced by 
the instability of the Kolmogorov flow 0911 . Marangoni instability of a two-layer system with a deformable interface 
ll^ . nonlinear development of zigzag instability of convection rolls EH, and even coarsening of ordered domains in 
oscillatory patterns governed by the complex Swift-Hohenberg equation 11421] . which describes oscillations in lasers 
El, El and optical parametric oscillations - El- 

For sake of simplicity, apply the boundary conditions (px - 4>xxx -Os.t x - ±1, then 


j(±l) 




(pdx — 0 


(the total length of domains of each phase is conserved). The Lyapunov functional (O decreases with time: 


F,^ 


f 


{(pxx + (p- (jy’fdx < 0. 


Using the boundary conditions, one can present equation (ISTT i in the form 


+ (pxx + (p - + h{t) = 0, d/(p,(x, t) = 


u: 


I-* - y\(litiy, t)dy. 


(32) 


(33) 


(34) 


Note that despite the energetical equivalence of phases (p - +1, an efficient field h(t) appears, which must to be 
determined self-consistently. The motion of a kink and an antikink towards each other would change the total lengths 
of domains of different phases, and hence it is not possible. Two kinks can move simultaneously towards an antikink 
placed between them, or kink-antikink pairs can move as a whole. The correlated motion of n domain walls with the 
centers at f = 1,..., n, sufficiently far from each other, is governed by the system 02.511 


= 16^e«l^sign(fy-^,) + 2(-l)m f- 


(35) 


2=1 


i*i 


supplemented by the conservation law 

n 

=0- (36) 

i=l 

For a kink-antikink pair, the attraction is compensated by the field h - 8 exp[-(^2 - V2], hence the domain walls 

are motionless. For a symmetric kink-antikink-kink triplet with coordinates of domain walls - -lit), ^2-0 and 
^3 = /(f), one obtains h — 0,1' - -8 exp(-/ V2), hence the annihilation time depends on /q = 1(0) as 

fo ~ -^loe'o^^ 

8 V2 


(cf. (HI). 
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2.2.2. Non-potential systems 

As an example of a non-potential system with a conservation law let us consider the convective Cahn-Hilliard 
equation, 

, D j 

^ ^ ^ ™ ( 37 ) 


+ {^xx + (l>-(f>)xx- —(4> )x^0, -oo <x<O 0 , 

which has been suggested to describe several physical processes, namely spinodal decomposition of phase separating 
systems in an external field ifs^ - 1^ . step instability on a crystal surface faceting of growing, thermody¬ 
namically unstable surfaces 15711 - 16ill , evolving nanofoams 16211 as well as dewetting of a thin film flowing down 
an inclined plane 16311 . That equation provides “a bridge” between the Cahn-Hilliard equation ([1]) (D = 0) and the 
Kuramoto-Sivashinsky equation (O (0 —> -20/D, D » 1). 

Stationary patterns 0 = 0(x) are described by the problem 


D(p^ DA 

-I- (0 — 0 )- — — - —oo < X < oo; A>0; X ^ ±oo. 


< oo. 


(38) 


For any D + 0, the set of solutions of equation (l38l l is incomparably more complex than that of the usual Cahn-Hilliard 
equation. One can easily And some exact solutions of the problem. The constant solutions. 


= 0± = + Va, 


(39) 


correspond to two stable phases. For domain walls, there exist exact solutions 05311 . one for a kink with A = A+ = 

1 -H D/ V2, _ 

0 = 0+(x) = 0+ tanh ■%:(x - ^), 0° = Jl + D/ V2, ^ = const, (40) 

V2 

and the other for an antikink with A - A- - \ - Dj V2, D < yfl. 


0_(x) = -0° tanh ^—{x - f), 0*! = Jl - D/ V2, f = const. 
V2 


(41) 


However, the set of stationary solutions is much more rich. Specifically, solution (l40l i is just one representative of a 
family of kinks 0+(x; A). In addition to the monotonic antikink (HTt . there exists also a discrete set of non-monotonic 
antikink solutions llh^l (that phenomenon is typical for models containing higher-order spatial derivatives, see Jh^]- 
ll^ l. The kink-antikink pair is formed by antikink (HTb and a representative of the fa mily of kinks with A = A_. 

If 0 < D < Do = V 2 / 3 . In that region, the coarsening is observed |55|, if^ . ifs^ . 16111 . Because of the asymmetry 
between kinks and antikinks, a kink-antikink pair moves spontaneously with a definite velocity V 2 (D, L). The most 
typical process observed by coarsening is the annihilation of domain wall triplets, when two kinks of the same sign 
approach with velocities +V 3 (D, L) the kink of the opposite sign situated between them. Exact expressions for V 2 (D, L) 
and V 3 (D, L) can be found in ll^ . In the limit of small D ll^ . 

V 2 (D, L) ~ V 3 (D, L) ~ -(D^ V 2 / 4 ) exp(-DL/2). 


Therefore, the coarsening law is logarithmic. 

For D > Do, the domain walls have oscillatory tails. That case will be discussed in the next section. 


3. Factors hindering coarsening 

In the present section, we discuss some typical situations where the system cannot reach a uniform state or another 
energetically preferred state by coarsening. 
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3.1. External inhomogeneities 

The motion of domain walls leading to annihilation can be stopped by inhomogeneity of the medium. Recall that 
we consider the phenomena in the absence of noise. Coarsening in inhomogeneous systems in the presence of thermal 
fluctuations is considered in . 

In a potential system, the domain wall would “prefer” the location where its energy will be smaller than in other 
locations. For example, let us consider the following modification of the one-dimensional Allen-Cahn equation iEbII : 

(p, = (pxx + [I + ef(x)](f> - (f>^. (42) 

At the leading order in e, the equation of motion for a domain wall of any kind is 

^ where Vih{^ = - I f(^ + y)cos\E^{y)dy. (43) 

3 dt 2 J_^ 

Specihcally, if the inhomogeneity has a (5-like shape, 

f(x) = -2Vod(x - X.), 


the interaction potential is 

V,7,(xo) = -Vocosh^^ 

V2 

Generally, the shape of the potential is a linear transformation of the inhomogeneity shape, according to (l43T l. 

If there are many domain walls and many inhomogeneities, the motion of domain walls is determined by the 
system of equations (fT4l i with the potential 

t/ = J]iT(^,-^,-i) + ^y,-,(x,), 

i i 

where W(^i - ^,_i) is determined by equation (fTsT l. and Vihixi) corresponds to (l43l l. Thus, the problem of hnding 
stationary solutions of (l42l i is equivalent to hnding equilibrium conhgurations of a chain of particles in the external 
potential 0^, interacting according to the law (fTsT l. This model resembles the well-known Frenkel-Kontorova model 
(see e.g. 1701]). 

As an example let us consider two domain walls with coordinates and ^2 which are near the distant d-shaped 
attracting inhomogeneities: 


fix) - -2yod(x - xi.) - 2yod(x - X 2 ,), Vo > 0, X 2 , > xi,. 

X 2 * - xi, - k » 1, 1^1 - xi,| = 0(1), 1^2 - X 2 *\ — 0(1). The equation of motion for the left domain wall is: 

-3 ~ 


2 V 2 - i 5 g (ft ft)V2_ V2yosinh^i—^ cosh 
3 dt a/2 


V2 


(44) 


The hrst term in the right-hand side of equation can be estimated as 16 exp(-/* V2). The minimum of the second term 
in the right-hand side of the equation is equal to -2 s/2Vol3 VS. Thus, we come to the conclusion that if 


2V2 
3 V3 


Vo > 16e-''^, 


(45) 


the domain wall will not be able to escape from the potential well created by the inhomogeneity. Hence, the coarsening 
will be stopped when the distances between the neighbor domain walls satisfy the inequality (l45T l. 

Similarly, in the case of domain walls pushed by the asymmetry of phases (see dSll, we find the criterion of 
pinning; 

3 V3, 


Vo > 


V 2 ' 


-h. 


If f(x) is a periodic function, the sequence of pinning sites (minima of the potential) filled by pinned domain walls 
can be regular (“commensurate patterns”) or irregular (“spatial chaos”) iIt^I . 
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3.2. Oscillatory tails of domain walls and stability of stationary patterns 

In the examples considered in Sec. 2, the domain walls are described by monotonic functions like (fTOl i. The mono¬ 
tonicity of the asymptotic behavior of the domain wall solution on the infinity leads to a sign-preserving (attracting) 
interaction between domain walls. Oscillatory tails of domain walls create a sign-alternating interaction potential. 
The domain walls can be captured near the potential minima, therefore stable patterns are formed due to pinning of a 
domain wall by an inhomogeneity created by another domain wall. 

As an example, let us consider the stability of periodic solutions of the convective Cahn-Hilliard equation (ITtT i. 
which satisfy the condition fix + 1) - fix). At large I, these solutions resemble periodic sequences of domain walls. 
Define the pattern wavenumber K - Injl. The normal disturbances of a periodic solution have the shape of a Floquet- 
Bloch function, f{x, t) - f{x) exp{ikx + crt), where f(x + l) - f{x), and k is a quasi-wavenumber, \k\ < Kjl. A periodic 
solution is always neutrally stable (cr = 0) with respect to a spatial shift, f{x) - fxix), k - 0. Therefore, a special 
attention has to be payed to potentially unstable longwave disturbances with small k. Their growth rate cr{k', K) can 
be presented as 

cr(k; K) - a-\{K)k + cr 2 (K)k^ H-... 

One can show that the sign of cr^{K) depends on the dependence of the squared pattern amplitude 


A = 


- r 

I Jo 


(p^(x)dx 


(see (l38T l) on the pattern wavenumber K = Itt/ I ifl^ . illtl, iS- If dAjdK < 0 for any K, which takes place for 
D < Do - V2/3, then crj(A') > 0 for any K, therefore all periodic solutions are unstable. That is compatible with 
the attractive interaction between domain walls. For D > Dq, the function A(/r) is not monotonic, and the extrema of 
AiK) separate the regions of a monotonic growth of longwave disturbances, cr\{K) > 0, and those of the oscillatory 


response of patterns to dilations and compressions, (PyiK) < 0 0741] . The region of oscillatory response can contain 
a subinterval of stable patterns (where cr 2 {K) < 0) 1^ . j^ . That is possible because of the alternating sign of 
the interaction between domain walls. On the boundaries of the stability interval, the patterns become unstable with 
respect to either longwave (phase) disturbances (see |75], or shortwave disturbances with k - Kjl, leading to a 
spatial period doubling. 

Because the potential of the interaction between domain walls has multiple minima, the distance between domain 
walls is not selected in a unique way. The pattern includes elements with “short” and “long” distances between the 
maxima that alternate in a rather irregular way S. 


3.3. Nonlocal interaction 

A specific kind of an arrested coarsening process has been found for equation ^ ifl^ . which can be written also 
as 

+ T4)) + fxx + f-f^+ h(t) = 0. (46) 

For small F the lowest density of the Lyapunov functional is achieved for patterns with wavelength Aopt = 6>(F '^^). 
According to the linear stability theory, the disturbance with largest growth rate has a wavelength Ac - 0{V). There¬ 
fore, one could expect that the energetically preferable, long-wave pattern will be developed from the initial short¬ 
wave pattern by coarsening. The coarsening process takes place indeed, but it is stopped when the wavelength reaches 
a much smaller value, /!„,■„ = 6>(ln(l/F)). The criterion of the pattern stabilization is similar to (l45T l. but now the 
stabilizing factor is the nonlocal interaction which is proportional to F. 


3.4. Pattern-induced pinning of a domain wall. 

The stability regions of a periodic pattern and a uniform state may overlap. In that case, the behavior of a domain 
wall between the pattern and the uniform state is crucial. Near the threshold of the instability creating short-wave 
patterns, where the width of a domain wall is large compared to the characteristic pattern wavelength, one can describe 
the dynamics of a domain wall using the envelope function approach if^ . if^ . In the framework of that approach one 
comes to the conclusion that a domain wall between the pattern and the uniform state moves with a constant velocity, 
which is proportional to the difference between the Lyapunov functional densities of the phases 114ill . However, the 
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influence of the underlying periodic pattern leads to some qualitative changes of the domain wall dynamics. First, 
the motion of the domain wall is an oscillatory process; during one period, one stripe is created or melted S, s. 
Secondly, because of the pinning effect, there is a hnite interval of the parameter value where the domain wall is 
motionless, i.e., a pattern and a uniform state coexist. Near the threshold of the pattern appearance, that interval is 
transcendentally small |8i], ED]. 

As an example, let us consider the competition and coexistence between patterns and uniform states for a system 
governed by the Swift-Hohenberg equation. 




dx^ 


+ 1 


which corresponds to the Lyapunov functional 

/{- 


FW 


IjP- + + 


1 


dx^ 


■ dx. 


(47) 


(48) 


That model was suggested for studying hydrodynamics fluctuations near the instability threshold ll82|] and used for 


modelling Benard convection 18311 . At e > 0, periodic patterns exist with wavenumbers k in the interval I - yfe <k^ 


1 + yfe, and they are stable in a subinterval A:_(e) < k < k+{e). At e > 1, constant nonzero solutions (p± — ± xfe- 1 
appear. At e > 3/2 they become stable with respect to small disturbances; kinks with oscillatory tails connect both 
stable uniform phases (f>± 08411 . 


The value of the Lyapunov functional density for the regular pattern with an optimal wavenumber is lower than 
that of the uniform state when e < em ~ 6.3 |84], S. Nevertheless, the domain wall between both states is immobile 
for much smaller values of e, e > ec ~ 1.7574. The reason is the self-induced pinning caused by the oscillatory 
asymptotic perturbation of the uniform state. Similarly, the pinning effect prevents the replacement of a pattern by 
a uniform state at e > The stability interval for a finite fragment of patterns sandwiched between semi-infinite 
regions of a uniform state slightly depends on the length of that fragment il% . Note that the noise activates the 
transition for a metastable state to a truly stable, energetically preferred, state 17911 . 

The coexistence of patterns and uniform states has been revealed for many pattern-forming systems (for a review, 
see isstl). 


4. Coarsening in two and three dimensions: curvature effects 

4.1. Phase separation 

First, let us consider domain coarsening for spatially uniform states in a system without the conservation law. 

In a two-dimensional (three-dimensional) potential system, the Lyapunov functional can be diminished without 
annihilation of a domain wall, just by diminishing its length (area). Let us consider a two-dimensional Allen-Cahn 
equation. 


<Pt - <Pxx + ^'yy + (p — (p 


(49) 


and present the isoline <p{x, y,t) — 0 (“a front”), which describes the center of a curved domain walls between the 
stable uniform phases, in the form y - H{x, t): <p{x,y, f) < 0 as y < H{x, t), (p{x,y, f) > 0 as y > H{x, t). One can 
show that the motion of the domain wall is determined, in the limit of small h and a small curvature of the front, by 
the equation [86], 1261. 

H, 3V2, 3V2, 

2 


Vtt:^ (i+i/2)3/2 


-h, QXV — K- 


-h. 


(50) 


where v is the normal velocity, and k is the curvature of the front. Specifically, in the case h 
same free energy), we get just the relation v = k, which is called curvature flow. 

As an example, let us consider a round droplet of the phase (p^ in the infinite sea of the phase (p 
in the case h - 0 the droplet radius is changed according to the law 


0 (both phases have the 
Because k - \ jR, 


R^{t) = R(0f ■ 
10 


2t. 
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The droplet collapses during the finite time f* « Rq/'^- The obtained life time of the droplet shows that the charac¬ 
teristic coarsening scaling is / ~ which is significantly faster than in the one-dimensional case. The same 

coarsening law is obtained in the 3D case. 

Moreover, even for the fractional Allen-Cahn equation. 


4), = + +h, 1 < y < 2, 


(51) 


the front motion is determined by a formula similar to (l50l l (up to numerical coefficients that depend on y), and the 
scaling I ~ is established on the late stage of coarsening (at the initial state, / ~ (9(f'^^), according to the 

scaling properties of the linearized equation) il. 

A more significant change of the front dynamics is produced by memory, when the temporal evolution of the order 
parameter is governed by the equation 

f' dF 

(p, - - I a{t - s) — (x,s)ds. (52) 

Jo 

Equation (l50l l is replaced by 

vr . 3V2, 


1 — av^ 


+ y V 


-hil-avy 


(53) 


where the constants a and y are determined by the Laplace transform of the kernel a(t - s) Il87n . 

In the presence of a conservation law, using the Cahn-Hilliard equation, one can find that the evaporation of a 
single round droplet of the phase in the infinite sea of the phase 0+ is governed by the equation |@1 


yit) = R^Q) - -o-t. 


where cr is a parameter corresponding to the effective surface tension of the domain wall. Hence, the scaling law 
I ~ Oyy is predicted. In the case of droplets of different sizes, the main mechanism of coarsening is the growth 
of big droplets (with a smaller curvature) at the expense of small droplets (with a larger curvature), which leads to 
“flattening” of the interphase boundary and hence the decrease of the Lyapunov functional (“Ostwald ripening”). 
Lifshitz and Slyozov and Wagner 0] have developed a kinetic theory of the phase separation in the limit of 
small concentration of the minority phase. A detailed description of that theory and its extensions can be found 
in 0211] . Here we mention some basic results. There exists a critical radius Rc(t) ~ such that smaller droplets 
evaporate by diffusion, while larger droplets grow by absorbing the matter through the majority phase. They have 
obtained a self-similar droplet radius distribution and found the law R{f) ~ for the characteristic droplet radius. 
Later, the latter law was confirmed, theoretically and numerically, for arbitrary concentrations of phases 191 l- j^ . 
Note that at shorter time after the beginning of the phase separation process, a scaling law R{f) ~ has been 
predicted and observed 1^ -1971. The same scaling laws are observed in the framework of a more general model. 


[M(0)(vV + 4 >- 0 )] 


(54) 


with a non-constant mobility function M{(f>). Lor instance, the crossover from R{t) ~ to R{t) ~ has been 
observed for M((^) = 1 - (p^ 1981. 

A nontrivial kind of Cahn-Hilliard equation has been derived for the description of flows in thin liquid films in 
the presence of disjoining pressure iIt^I - IdiIi . Here the conservation law is the the conservation of the liquid volume, 
while “the phases” are a macroscopic film and a mesoscopic “precursor film”. The Cahn-Hilliard equation describes 
the decomposition of a film into droplets connected by the thin precursor film. The coarsening of droplets is due 
to the growth of large droplets at the expense of small ones and because of the motion of droplets leading to their 
coalescence. In the framework of the standard model ll49n . the coarsening law is N(t) ~ where N{t) is the 

number of droplets. Generalizations of that model leading to different coarsening rates can be found in 0211. 

_ A specific kind of coarsening takes place if there are more than two thermodynamically equivalent phase s 

EM, e.g., because of different possible orientations of the spin (the list of examples can be found in jlOOll l. The 
problem can be modelled by means of an overdamped sine-Gordon equation similar the Allen-Cahn equation O but 
with the potential V{(f>) — cos - 1; the stable equilibrium phases correspond to the potential maxima, (p — 2nmlp, 
m - 0,1,..., p - 1. Lor p - 2, the coarsening is similar to that for the Allen-Cahn equation in any dimension. If 
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p > 2, the coarsening in ID is determined by the exponentially weak interaction of domain w alls, which can be now 
either attractive or repulsive. For p > d + 1, a logarithmic rate of coarsening is also predict ed [llOOl] . 

Coarsening in non-potential systems was studied using isotropic ll^ and anisotropic i59ll generalizations of the 
convective Cahn-Hilliard equation. 

A numeri cal an alysis of coarsening versus pattern formation in non-potential pattern-forming systems has been 
carried out in Cml. 


4.2. Pattern ordering 

The naturally appearing patterns usually contain numerous defects 01021] . Specifically, the pattern may have a 
multidomain structure. One can distinguish between two kinds of domain walls in periodic patterns irill . The first 
kind of domain walls separates patterns of different symmetry, which are generally not energetically equivalent (e.g., 
stripe patterns and hexagonal patterns). In that case, the domain wall tends to move expanding the energetically 
prefered domain, but it can be stopped by pinning, as it was explained in the previous sections 10. The second 
kind of domain walls separates patterns of the same type but with different orientations or different values of the 
wave length. Domains of different orientations appear spontaneously or are created by the side walls in a finite region 
1 103 1. Different scenarios of domain wall evo lutio n are possible 0: (i) the domain wall can be a source of a 
waven umbe r selection, similarly to a side wall llOdll or a ramp smoothly matching pattern region with a subcritical 
region illOSll : (ii) the domain wall can be destroyed by an intri n sic in stability; (iii) it can spread and smoo th dow n. 

Besides d omai n walls , patterns contain dislocations llOhll . 1107], coupled pairs of dislocations |108L j 109L and 
disclinations Q-ilil. Their motion is also a significant factor of the ordering in periodic patterns lld^l . 

While some specific phenomena related to the dynamics of defects in patterns have been a subject of a theoret¬ 
ical analysis, the full picture of pattern ordering is studied mostly by means of numerical simulations of dynamical 
equations (possibly with noise) or experimentally. 

In isotropic systems with symmetry 0^-0 (e.g., for (|2]) or (l47l) '). stripes of different orientations are generated. 
That allows a simplified descr iptio n of patterns by the phase field <p(x,t), <p(x,t) - 0o cos[i/'(x, f)], with the local 
wavevector k(x, t) - Vi/'(x, t) 1 11 311 . The numerical simulations reveal power laws for the grow th of domains and 
eliminati on of domain walls, dislocations and disclinations for the original equation (|47]) ll 1411 and for the phase 
equation ill 1-511 . Computations carried out for (l2]i show that on the background of the final state, which is a unidomain 
structure, the orientational two-time correlation function has a power-law asymptotics, while the spatial two-point 
correlation function is subject to a transition from a power law to an expo nential law with time Eil. 

If the symmetry 0 —> -0 is violated (due to an external field ll ITH or a cubic term in the free energy density 


1 1 18l1 I. a competition between stripes and hexagons takes place. As an example, let us mention ordering in patterns 
governed by a generalized Swift-Hohenberg equation 


01 


= = J [- 2 * *?* 


+ ^0^ + ^[(V^ -H l)0]^^^fx, X = (x,y). 


(55) 


Changing s, on e can arrange a transition between stripes and hexagons and vice versa. The analysis of the transition 
has been done im by studying the structure factor S (k, t) - (|0(k, f)P), where 0(k, t) is the Fourier transform of 
the order parameter (pix, t), and () denotes ensemble averaging. One has found different scaling laws for ordering the 
stripes, for the growth of hexagonal domains due to the stripe-to-hexagon transition, and for the growth of stripes from 
a disordered hexagonal patterns. 

Orientational ordering in hexagonal patterns has been studied experimentally in im and numerically in il20ll . 
using the modification of equation (l2]i with broken inversion theory. Note that the problem of t he or ientational ordering 
is related to the problem of coarsening in the system with degenerate phases studied in [99], lIlOOll . 

In conclusion, we have reviewed basic models and effects related to coarsening in pattern forming systems. Re- 
centl y, investigations of more complex cases have b een i nitiated, e.g., domain coar senin g in an oscillatory patterns 
II 121 II . pattern coarsening in time-dependent domains 012211 . and patterns in networks 11230 . These subjects are beyond 
the scope of the present review. 
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